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95°F-days  to  783°F  days.  Surface  freezing  indexes  in  asphaltic  concrete  pave- 
ments ranged  from  59°F-days  to  400°F-dgys,  and  for  Portland  cement  concrete 
pavements  ranged  from  57  F-days  to  388  F-days.  All  freezing  indexes  were  con- 
siderably greater  during  the  1976-1977  winter.  'Measured  maximum  frost  penetra- 
tion depths  were  compared  with  computed  maximum  depths  for  each  of  30  test  sites 
The  modified  Berggren  equation  was  used  to  compute  the  maximum  frost  depth  at 
each  site.  Measured  maximum  frost  depths  ranged  from  20.5  in.  to  52.0  in., 
while  computed  maximum  values  ranged  from  14.0  in.  to  61.0  in.  The  mean  dif- 
ference between  observed  and  computed  maximum  frost  penetration  depths  was 
3.8  in.  ‘Maximum  frost  penetration  depths  were  computed  for  hypothetical  pave- 
ments with  open-graded  drainage  at  four  of  the  test  sites.  It  was  concluded 
that  open-graded  drainage  layers  would  not  significantly  change  the  frost  pene- 
tration beneath  highway  pavements  in  New  Jersey.  It  was  recommended  that  test 
pavements  be  installed  to  verify  the  computations. 
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CONVERSION  FACTORS:  U.S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENT 

These  conversion  factors  include  all  the  significant  digits 
given  in  the  conversion  tables  in  the  ASTM  Metric  Practice  Guide 
(E  380),  which  has  been  approved  for  use  by  the  Department  of  Defense. 
Converted  values  should  be  rounded  to  have  the  same  precision  as  the 
original  (see  E 380). 
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IMPROVED  DRAINAGE  AND  FROST  ACTION  CRITERIA 
FOR  NEW  JERSEY  PAVEMENT  DESIGN  PHASE  II 
(DATA  ANALYSIS) 

Richard  L.  Berg 
INTRODUCTION 

In  early  1975,  engineers  from  the  New  Jersey  Department  of  Trans- 
portation (NJDOT)  visited  CRREL  to  discuss  a program  they  were  embarking 
upon  with  assistance  from  the  Federal  Highway  Administration.  NJDOT 
was  interested  in  including  open-graded  drainage  layers  in  their  pave- 
ment, but  before  doing  so  chose  to  determine  analytically  the  effect  of 
frost  penetration  beneath  pavements  when  using  an  open-graded  drainage 
layer.  In  early  July  1975  NJDOT  issued  a contract  to  CRREL  to  conduct 
the  study,  the  principal  goal  of  which  was  to  analytically  determine  the 
impact  of  an  open-graded  drainage  layer  on  the  depth  of  frost  penetra- 
tion beneath  pavements  in  New  Jersey. 

The  final  report  on  that  study  (Berg  and  McGaw  1978)  was  submitted 
to  NJDOT  In  June  1978.  In  the  previous  study,  the  modified  Berggren 
equation  was  used  to  estimate  the  maximum  seasonal  frost  penetration 
depths  beneath  several  pavement  profiles  In  New  Jersey.  Mean  and  design 
freezing  indexes  were  developed  for  Newton,  Trenton  and  Altantic  City, 

New  Jersey,  and  used  in  the  calculations.  Surface  transfer  coefficients, 
i.e.,  n-factors,  were  determined  from  air  and  surface  temperature  data 
supplied  by  NJDOT.  They  also  furnished  aggregates  for  thermal  conduc- 
tivitv  tests  on  base  and  subbase  course  materials  and  paving  materials. 
Aggregates  for  the  proposed  stabilized  and  unstabilized  drainage  layers 
were  also  Included  for  thermal  conductivity  tests. 

The  current  study  was  initiated  in  December  1976  and  was  to  extend 
through  September  1977.  In  October  1977  the  contract  was  extended  to 
30  March  1978. 

The  study  had  the  following  objectives: 

a.  To  determine  the  probable  magnitude  of  frost  penetration 
depths  at  specific  sites  in  New  Jersey. 

b.  To  compare  calculated  frost  penetration  depths  with  measured 
depths  to  verify  the  prediction  capability  of  the  modified  Berggren 
equation  in  New  Jersey,  and 

c.  To  determine  whether  there  may  be  special  or  unanticipated 
adverse  affects  due  to  application  of  open-graded  drainage  layers. 


THERMAL  ANALYSIS 


In  the  previous  studv  (Berg  and  McGaw  1978),  pages  23-34,  the 
modified  Berggren  equation  was  used  to  calculate  maximum  frost  pene- 
tration depths  in  selected  New  Jersey  pavement  profiles.  At  the  sugges- 
tion of  CRREL  engineers,  NJDOT  engineers  and  technicians  installed  frost 
tubes  for  monitoring  frost  penetration  depths  and  access  tubes  for 
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monitoring  subsurface  moisture  conditions  with  a nuclear  probe.  CRREI.'s 
N.IDOT  colleagues  collected  and  tabulated  data  from  the  test  sites  and 
sent  the  data  to  CRREL.  The  remainder  of  this  report  describes  infor- 
mation from  the  field  test  sites  and  gives  calculations  of  frost  penetra- 
tions at  the  sites. 


Field  Test  Sites 


During  the  1976  summer  and  fall  N.IDOT  installed  frost  tubes  and 
access  tubes  for  a nuclear  moisture  sensor  at  30  locations  in  the 
northern,  central  and  southern  portions  of  New  Jersey.  Figure  la  shows 
the  locations  of  the  test  sites.  During  the  1976-1977  winter,  N.IDOT 
personnel  measured  frost  penetration  depths  and  moisture  contents  at 
these  sites  two  to  three  times  per  week.  The  frost  penetration  and 
moisture  content  data  were  tabulated  bv  N.IDOT  personnel  and  mailed  to 
CRREL.  At  CRREL  the  data  were  analyzed  and  used  in  the  computations  of 
frost  penetration  depth.  All  of  the  frost  penetration  data  were  plotted 
as  were  measured  moisture  contents  at  selected  depths;  these  are  presented 
in  Appendix  A. 

Figure  lb  shows  three  locations  in  the  state  (Rockaway,  Bedminster 
and  Bordentown)  where  air  and  pavement  surface  temperatures  were  measured 
simultaneously.  Air  and  surface  temperatures  were  measured  at  the  same 
three  sites  during  the  1975-1976  winter  and  their  respective  freezing 
indexes  were  reported  previously  (Berg  and  McCaw  1978),  pages  5-6. 

Daily  maximum  and  minimum  air  and  pavement  surface  temperatures  were 
tabulated  bv  N.IDOT  personnel  and  sent  to  CRREL  where  freezing  indexes 
and  n-f actors  were  determined.  At  each  of  the  three  sites,  surface 
temperatures  were  measured  on  both  asphaltic  concrete  and  Portland 
cement  concrete  pavements. 

Six  sites  were  choosen  bv  NJDOT  (1977)  for  detailed  thermal  analysis 
using  the  modified  Berggren  equation.  At  these  sites,  different  para- 
meters in  the  modified  Berggren  equation  were  varied  to  evaluate  their 
influence  on  computed  frost  depths.  After  completion  of  these  calcula- 
tions, frost  penetration  depths  were  computed  and  compared  to  measured 
values  for  the  remaining  24  sites.  The  six  sites  chosen  tor  detailed 
analysis  were  Site  Numbers  6,  9,  11,  18,  22  and  29  (Figure  la).  The 
thickness  and  properties  of  various  materials  in  the  pavement  protiles 
are  discussed  subsequently. 

Surface  Transfer  Coefficients 

Berg  and  McCaw  (1978)  discussed  surface  transfer  coefficients,  n- 
factors,  which  have  been  observed  in  various  geographical  locations  in 
the  United  States,  Canada  and  Norway.  As  previously  stated  the  New 
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Jersey  Department  of  Transportation  measured  surface  and  air  tempera- 
tures during  the  1975-1976  and  1976-1977  winters. 


Data  for  the  1975-1976  winter  were  reported  previously  (Berg  and 
McCaw  1978),  and  freezing  indexes  and  n-factors  measured  during  the 
1976-1977  winter  are  shown  in  '•''able  1.  At  all  three  sites  the  air 
freezing  index  and  length  of  the  air  freezing  season  were  greater  in  the 
1976-77  winter  than  in  the  previous  winter.  The  same  is  true  for  the 
Portland  cement  and  asphaltic  concrete  pavement  surface  freezing  indexes. 


Table  1.  Freezing  indexes  and  n-factors 
1976-1977  winter. 


Air 

Portland 

cement 

concrete 

Asphal 

tic  concrete 

Season 

Index 

Season 

Index 

n- 

Season 

Index 

n- 

Site 

days 

°F-days 

days 

°F-days 

factor 

days 

°F-days 

f actor 

Bordentown 

46 

316 

41 

93 

0.29 

44 

181 

0.57 

Bedminster 

51 

446 

52 

388 

0.87 

52 

400 

0.90 

Rockaway 

84 

783 

72 

295 

0.  38 

70 

304 

0.39 

Table  2 summarizes  n-factors  for  the  asphalt  concrete  and  Portland 
Cement  concrete  pavements  during  1975-76  and  1976-77  winters.  For 
Bordentown  and  Rockaway  the  n-factors  were  generally  slightly  lower 
during  the  1976-1977  winter  than  observed  in  the  1975-1976  winter.  The 
exception  was  the  n-factor  for  the  Portland  cement  concrete  test  section 
at  Bordentown  which  was  considerably  less  during  the  1976-1977  winter 
than  observed  during  the  1975-1976  winter.  Data  from  Bedminster  exhibited 
a different  behavior.  During  both  winters  the  n-factor  for  asphaltic 
concrete  was  considerably  higher  at  Bedminster  than  observed  at  the 
other  two  sites.  Data  from  the  Bedminster  site  also  indicated  that  the 
n-factor  for  the  asphaltic  concrete  pavement  was  slightly  higher  during 
the  1976-1977  winter  than  measured  during  the  1975-1976  winter.  Although 
the  n-factor  over  the  Portland  cement  concrete  pavement  at  Bedminster 
agreed  quite  well  with  the  observations  at  the  other  two  sites  during 
the  1975-1976  winter,  it  was  considerably  higher  during  the  1976-1977 
winter . 

Table  2.  n-Factors  from  New  Jersey  sites. 


Asphaltic  concrete  Portland  cement  concrete 


Site 

1975-1976 

1976-1977 

1975-1976 

1976-1977 

Bordentown 

0.62 

0.57 

0.60 

0.29 

Bedminster 

0.84 

0.90 

0.54 

0.87 

Rockaway 

0.46 

0.39 

0.40 

0.38 

Use*  0.55 

Use*  0.45 

*Value  used  in  computation  of  maximum  frost  depth  using  the  modified 
Berggren  equation  if  the  surface  freezing  index  was  not  measured  at 
a nearby  location. 


Air  and  surface  temperature  data  information  from  Bedminster  show 
no  anomalous  behavior  during  either  the  1975-1976  or  1976-1977  winters. 

As  part  of  the  current  study,  however,  air  freezing  indexes  during  the 
1976-1977  freezing  season  were  computed  for  several  other  locations  in 
New  Jersey  (Table  3).  These  data  indicate  that  air  temperatures  n assured 
at  Bedminster  were  significantly  higher  than  those  meas'red  in  nearby 
communities,  thereby  cr  sing  a lower  freezing  index  at  Bedminster.  Air 
freezing  indexes  at  Long  Valley  and  Somerville  averaged  795  degree  days 
for  the  1976-1977  winter,  whereas  the  air  freezing  index  measured  at 
Bedminster  was  only  446  degree  days. 

The  lengths  of  the  freezing  seasons  were  significantly  different  at 
the  three  sites  also. 


Table  3.  Weather  stations,  air  freezing  indexes 
and  lengths  of  freezing  season  used 


in 

the 

thermal 

analysis . 

Site 

No 

Weather 

station* 

Location** 

Air  freezing 
index, 
°F-days 

Length  o 
freezing  : 
son,  dayi 

6 

Trenton 

12 

mi  W 

446 

51 

9 

Bordentown 

at 

site 

316 

46 

11 

Long  Valley 

9 

mi  N 

800 

85 

11 

Somerville 

9 

mi  SE 

790 

85 

18 

Canoe  Brook 

5 

mi  SE 

771 

86 

18 

Morris  Plains 

4 

mi  NW 

764 

86 

22 

Boon ton 

1 

mi  S 

822 

86 

29 

Newton 

7 

mi  W 

1075 

86 

* National  Oceanic  and  Atmospheric  Adminstration  stations  except 
Bordentown  where  temperatures  were  obtained  by  NJDOT. 

**  Approximate  distance  and  direction  from  test  site. 


At  both  Long  Valley  and  Somerville  the  freezing  season  was  85  days, 
but  the  season  observed  at  Bedminster  was  only  51  days.  Using  the 
average  air  freezing  index  (795°F-days)  from  Long  Valley  and  Somerville 
with  the  surface  freezing  indexes  (388  and  400°F-days)  measured  at 
Bedminster,  the  calculated  n-factors  of  0.49  and  0.5°F  for  the  Bedminster 
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site  agree  much  more  closely  with  values  observed  at  the  other  two 
sites.  The  effect  of  variations  in  the  n-faotors  and  the  influence  of 
these  variations  were  discussed  previously  (Berg  and  McCaw  1978). 

Conditions  used  in  the  computations  for  each  of  the  sites  are 
summarized  subsequently,  and  the  n-faetor  used  when  the  surface  temper- 
ature was  not  measured  near  the  field  test  site  is  shown  in  Table  2. 


Thermal  Properties 

For  the  calculations  in  this  report,  thermal  properties  were  deter- 
mined by  methods  discussed  by  Berg  and  McCaw  (1978).  Moisture  contents 
for  each  of  the  field  test  sites  were  estimated  from  measurements  made 
bv  N.TDOT  personnel  and  furnished  to  CKRK1.  in  tabulated  torm.  The  moisture 
content  and  the  estimated  density  influence*  all  of  the  thermal  proport ios, 
i.e.,  thermal  conductivity,  volumetric  heat  capacity  and  volumetric 
latent  heat  of  fusion. 

The  volumetric  heat  capacity,  in  BtuA-u  it  °F,  is  determined  from 
the  following  equation: 

C = >d  ( + (0.75)  (w/100)|  (1) 

where  y = the  dry  unit  weight,  lb/cti  ft 

d 

c - the  specific  heat  of  dry  soil,  Btu/lb  °F 
w *=  the  moisture  content  , Z drv  weight 

The  volumetric  la t on t heat  ol  fusion,  in  Btu/cu  It,  is  determined 
from  equation  (2): 

I.  = 144  y , (w/100)  (2) 

(1 

The  thermal  conductivity  was  determined  by  the  procedure  outlined 
bv  Berg  and  McCaw  (1978): 

a.  For  the  unknown  soil,  determine  or  estimate  the  drv  density  and 
moisture  content  by  drv  weight. 

b.  Using  the*  equations  below  or  Figures  and  1 find  (version's 

value  of  thermal  conductivity,  K , lor  a soil  having  the  same  density 

s 

and  moisture  content. 


For  unfrozen  sandy  soils 


K - (0.7  log  <w)  t 0.411  10°'(1  \|  ( 1) 

and  lor  frozen  sandy  soils 

K - 0.07b  (10)°'°M\l  + 0.017  ( 10'»<)'n,4f’\l  (w)  (4) 

s 

where  K = thermal  conductivity,  Btu/so  tt  hr  °F/ln 

s 

> » drv  density  of  the  soil,  lh/eu  It 

u 

w “ moisture  content,  t drv  weight 


e.  Knter  the  appropriate  curve  tor  Figure  4 or  S with  the  per- 
centage passing  the  700  sieve  for  the  New  .lersov  soil  and  tlnd  the 
corresponding  ratio  of  K/K  . 

d.  Multiply  the  Kersten  value,  K , hv  the  ratio  determined  in  step 
c and  the  resulting  value  should  he  within  +107,  of  the  thermal  conduc- 
tivity for  the  New  Jersey  so i 1 . 

For  the  pavement  material  and  the  stabilized  and  unst ah i 1 i zed 
drainage  layers,  thermal  conductivity  values  used  in  tin'  calculations 
were  those  determined  from  thermal  conductivity  tests  conducted  during 
the  previous  study  (Berg  and  Me  Caw  |D78). 


Frost  Depth  Calculations 

Thicknesses  and  representative  values  I or  the  moisture  contents  ol 
the  various  paving  materials  were  obtained  from  data  provided  hv  the 
N.T00T.  Densities  for  each  of  the  materials  were  obtained  ! rom  data  used 
in  the  previous  study  (Berg  and  McOaw  10781,  and  t rom  N.1P0T  (1078). 

Densities  and  moisture  contents  were  used  with  the  methods  described 
above  to  estimate  the  thermal  proport ies  tor  each  si'll  liver;  results 
are  shown  in  Table  4. 

As  previously  stated  a detailed  analysis  ol  the  influence  ot  various 
parameters  was  performed  at  six  ol  the  selected  test  sites.  Table  S 
contains  initial  and  upper  boundary  condit ions  used  at  those  six  sites. 
Also  shown  In  Table  S are  the  frost  penetration  depths  computed  using 
each  set  of  initial  and  boundary  condit  ions  in  conjunct  ion  with  the  soil 
and  pavement  proport ies  shown  itt  Table  4.  The  t rest  dent  i was  calculated 
using  the  modified  Berggren  equation,  a complete  explain  ion  ot  which  is 
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(%)  Passing  No  200  Si«v«  (NJ  Motsnolt) 

Figure  4.  Comparison  of  thermal  conductivity  (K)  of  New  Jersey  soils 
with  Kersten's  values  (K^)  for  sandy  soils,  unfrozen. 


(%)  Passing  No  200  Si»*e  (NJ  MoUnols) 


Figure  S.  Comparison  ot  thermal  conductivity  (K1)  ot  Now  Jerscv  soils 
with  Kersten's  values  (K  1 for  sandv  soils,  frozen. 


Btu  . . Btu  . . Btu  . Btu  . , Btu  . , Btu 


S-M-G  130.  11.4  1.0  2.30  1.38  33.  2134.  1760.  1385 

M-c  130.  13.8  2.0  2.40  1.44  27.  2583.  1928.  1273 
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C-S-M  125.  6.9  1.0  1.48  - 28.  1242. 

G-S-M  125.  4.7  3.0  1.20  - 26.  846. 


Symbols  are  for  the  following  materials: 
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Table  5.  Boundary  conditions,  initial  conditions  and  thermal  properties  varied  for 
sites  given  a detailed  thermal  analysis. 


Mean  Annual  Temperature 


given  on  pages  23-25  in  the  earlier  study  by  Berg  and  McCaw  (1978).  At 
sites  9 and  11  surface  temperatures  were  measured;  thus  surface  freezing 
indexes  were  computed  direct Iv  from  temperature  measurements.  At  the 
other  four  sites  surface  freezing  indexes  were  obtained  bv  multiplying 
tiie  air  freezing  index  by  the  applicable  n-factor  from  Table  2. 

Initially  frost  depth  computations  were  made  using  values  of  thermal 
conductivity  (K  ) obtained  from  Kersten's  (1949)  equations  and  mean 
annual  temperatures  (MAT)  which  were  thought  to  be  representative  of  the 
sites,  see  Solution  1 for  Sites  b,  11  and  29  (Table  5).  The  resulting 
computed  frost  depths  (13,  20  and  12  In.)  wore  much  lower  than  the 
observed  values  (38,  38  and  52  in.)  in  all  cases.  The  author  felt  that 
using  the  average  soil  temperature  existing  Immediately  prior  to  onset 
of  freezing  would  be  more  realistic  than  using  the  mean  annual  tempera- 
ture which  is  typically  used  in  modified  Berggren  equation  solutions. 
Subsurface  soil  temperatures  are  measured  at  unpaved  areas  near  New 
Brunswick  and  Newton  and  are  reported  in  the  monthly  weather  summaries 
(NOAA  1976). 

Data  from  these  sites  for  20  December  1976  indicated  that  tempera- 
tures to  a depth  of  40  in.  ranged  from  1 )°K  to  43°F.  As  snow  was  allowed 
to  accumulate  on  these  plots  and  as  the  surface  was  vegetated  (or  bare) 
rather  than  paved,  one  would  expect  that  the  amount  of  cooling  beneath 
the  paved  areas  would  probably  have  been  more  extensive  than  beneath  the 
observed  plots.  Therefore  calculations  were  made  using  an  initial  soil 
temperature  of  35°F.  For  the  three  sites  (Nos.  6,  11,  29)  the  only 
difference  between  Solution  l and  Solution  2,  Table  5,  is  the  change  in 
MAT.  The  computed  frost  depths  (28,  34  and  49  in.)  agreed  much  more 
closely  with  measured  depths  with  this  change.  It  was  felt  that  35°F 
was  a reasonable  assumption  for  an  average  temperature  within  the  sea- 
sonal frost  zone  prior  to  freezing.  This  value  then  was  used  tor  t hi’ 

MAT  in  all  of  the  other  calculations. 

At  Site  No.  6,  the  n-factor  was  increased  to  0.9  tor  Solution  ), 
and  the  resulting  calculated  frost  depth  (15  in.)  was  slightly  less  than 
the  measured  depth  (38  in.).  At  Site  9,  the  measured  surface  freezing 
index,  93°F— days,  was  used  in  Solut ion  1 . The  computed  I rest  depth  (14 
in.)  was  about  one-third  ot  the  measured  depth  (40  in.).  lor  Solutions 
2 and  3 at  Site  9,  the  air  freezing  Index  and  n-factors  for  PCC  and  AC 
pavements,  respectively,  were  used.  Resulting  computed  f rost  depths  (18 
and  20  in.)  were  onlv  about  50%  of  the  observed  depth.  Consequently 
additional  computer  trials  with  various  surface  freezing  indexes  were 
used  in  an  attempt  to  more  closely  approximate  the  observed  depth  of 
freezing.  Solutions  7—10  at  Site  9 contain  these  results  and  indicate 
that  a surface  freezing  index  of  about  tViO  °l'-days  would  have  been 
required  to  produce  the  observed  t rost  depth  using  these  soil  proport  ios. 


This  exercise  shows  that  a surface  index  in  excess  of  the  air  index 
(316°F-days)  would  have  been  necessary,  an  untenable  situation.  Solution 
11  indicates  that  if  the  latent  heats  of  fusion  are  further  reduced  to 
L?  (Table  4)  and  the  thermal  conduct ivites  are  reduced  to  reflect  more 
representative  values  (K)  anticipated  for  New  Jersey  soils,  a surface 
freezing  index  of  446  °F-days  would  result  in  41  in.  of  frost  penetra- 
tion. Using  with  Kersten's  thermal  conductivity  value  (K  ) and  a 
surface  freezing  index  of  316  °F-davs  resulted  in  a predicted  39  in.  of 
frost  penetration.  It  appears  that,  for  frost  to  penetrate  to  the 
measured  depth  of  40  in.,  a significant  amount  of  unfrozen  moisture  must 
exist  in  the  soils  at  Site  9.  Unfortunately  no  information  is  available 
concerning  the  amount  of  unfrozen  moisture  in  the  New  Jersey  soils; 
however,  the  existence  of  unfrozen  water  in  soils  at  temperatures  below 
32°F  is  a well  accepted  phenomenon.  A typical  relationship  between 
unfrozen  water  content  and  below  freezing  temperatures  is  shown  in  Table 
6,  from  Lovell  (1957).  The  Influence  of  unfrozen  water  content  is 
particularly  significant  in  the  finer  grained  soils  at  temperatures  Just 
below  the  freezing  point. 

Table  6.  Comparative  phase  composition  at 
selected  temperatures. 


Unfrozen  moisture  content  w^  (Z) 

26 ,6°F  5°F  - 1 3 ° F Ratio 


Clayey  Silt 

Silty  Clay 
Clay 

4.4 

11.05 

23.6 

3.1 

8.0 

15.8 

2.8 

7.2 

13.9 

1.55 

1.53 

1.70 

Percentage  moisture 

w-w 

frozen  ~ 100 

w 

26.6° 

F 5°F 

- 1 3 ° F 

„ , - 1 1 ° 1 ' 
Ratio  TbVfT’K 

Clayey  Silt 
Silty  Clay 
Clay 

73.4 

41.8 

16.3 

81.0 

58.0 

44.0 

82.9 

62.1 

50.6 

1.11 

1.48 

3.11 

as 


For  each  of  the  6 s 
follows : 


a) 

Solut ion 

4 

vised 

b) 

Solut ion 

5 

used 

ltes.  Solutions 

K and  L, 
s 1 

K and  Lj 


4, 


5 and  b were  identical  except 
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c)  Solution  h used  K ami  1., 

Frost  depths  computed' using  t liosu  three  conditions  did  not  differ 
groat  lv  except  at  Situ  IS.  Thu  pavement  at  tills  situ  was  const  mot  ed 
ovor  a low-molsturo  ooutont  slltv  satnl  subgrnde  aiul  tho  reduct  ion  In  t ho 
latent  heat  of  fusion  in  Solution  h oausoil  a significantly  larger  increase 
in  computed  frost  penetration  than  at  the  other  five  sites. 

At  Sites  11  ami  22.  the  effects  of  longer  or  shorter  freezing 
seasons  with  the  same  freezing  index  were  studied.  The  onlv  difference 
between  Solutions  2 and  1 at  Site  11  was  the  length  ot  the  tree/.  Ing 
season.  The  shorter  season  and  its  colder  average  temperatures  caused  a 
slight  Increase  in  t rest  penetration  ( 14  vs.  1.'  in.).  At  S i t u 

Solutions  1 and  ( were  identical  except  that  the  tree/,  lug  season  was 

shorter  in  Solution  ).  Computed  I tost  depths  are  the  same  ( '•»  in.) 
in  both  solutions.  From  study  ot  these  two  sites.  It  is  concluded  that 
t ho  length  ot  the  freezing  season  does  not  significantly  affect  computed 
trust  depths.  This  conclusion  mav  not  he  valid  tor  freezing  seasons 
shorter  than  about  40  davs,  however. 

A careful  study  ot  Solut Ions  4,  S and  h tor  each  s i t e and  conpat I 

son  ot  results  using  the  same  values  t or  the  ait  trou/ing  index,  n 

factor,  length  ot  freezing  season  and  mean  annual  tumputaturu  (MAT)  with 
thermal  properties  K and  1.  indicates  that  no  significant  increase  in 
accuracy  ot  the  computat ions  is  achieved  hv  modifying  the  thermal  con- 
ductivity or  the  latent  heat  ot  fusion.  lahlc  summarizes  this  com 
pari  son. 

The  mean  difference,  x.  was  computed  t tom: 


where  x4  is  the  difference  for  the  ith  computation  and  n is  the  numbei 
of  c output  at  ions.  The  standard  doviat  ion,  s,  listed  in  the  table  t s tin 
standard  deviation  from  the  mean  ditiorence  and  is  dot  etmf  nod  Irom: 


In  comparing  computed  and  measured  I rest  penetration  depths,  the  mean 
difference  and  the  standard  deviation  Irom  the  mean  using  K,  and  1(  were 
i.H  In.  and  + 1 1 . h in.  Using  K and  1 the  values  were  '.0  In.  anil  + 0 • 0 
in.,  and  using  K and  !.,the  values  were  1.7  in.  and  * in. 


i 


1 

Because  the  computations  using  K and  L are  much  more  direct  and 
make  use  of  the  Modified  Berggren  equation  without  alteration,  these 
values  were  used  in  the  computations  for  the  other  24  sites.  Initial 
and  upper  boundary  conditions  for  the  24  sites  are  shown  in  Table  8; 
computed  frost  depths  are  shown  in  Table  9.  The  mean  difference  between 
computed  and  measured  results  was  3.8  in.  with  a standard  deviation  of 
+7.8  in.  It  is  observed,  however,  that  the  mean  difference  for  the 
first  ten  sites  was  8.9  in.  with  a standard  deviation  of  +7.8  in.  (If 
Site  9 is  omitted  the  mean  error  and  standard  deviation  reduce  to  b.9 
in.  and  + 5.5  in.,  respectively.) 

For  the  remaining  18  sites  (11-30),  the  mean  difference  between 
measured  and  computed  frost  depths  was  1.0  in.  with  a standard  deviation 
ot  + t> . 2 (n.  The  first  ten  sites  have  surface  freezing  indexes  less 
than  about  250  °F-davs,  indicating  that  the  reliability  of  the  computed 
frost  depths  for  the  lower  freezing  Indexes  is  less  than  for  the  higher 
Indexes.  Since  the  first  10  sites  are  south  of  Somerville  (Figure  la), 
it  Is  recommended  that  a lower  reliability  be  placed  on  computed  frost 
depths  south  of  this  location. 

Computed  frost  penetrations  mav  have  agreed  slightly  better  with 
measurements  had  a combination  of  K and  I.  or  K and  I. ^ been  used;  how- 
ever. based  on  results  presented  in'sTable  /,  improvement  would  have  been 
relatively  minor.  More  accurate  computations  would  probably  have  been 
made  if  the  average  soil  temperature  at  the  beginning  ot  the  freezing 
season  were  known  and/or  it  the  surface  freezing  Index  were  measured  at 
each  site.  Calculated  frost  depths  (Table  5)  indicate  that  signit  leant 
changes  in  frost  depth  occur  with  relatively  small  changes  in  surface 
freezing  index. 

Moulton  and  Pubbe  (19b8),  in  a similar  studv,  used  a slightly  dit- 
ferent  form  of  the  modified  Berggren  equation,  but  their  results  were 
similar  to  those  presented  in  this  report.  Thev  concluded  that  accurate 
calculated  results  would  be  obtained  provided  reasonable  values  ot  the 
controlling  parameters  were  used.  They  state,  "These  parameters  include 
the  pavement  surface  freezing  index,  an  init  ial  soil  temperature  equal 
to  the  average  soil  temperature  at  the  start  ot  pavement  surtace  freezing 
and  reasonable  accurate  values  of  the  thermal  properties  ot  the  pavement 
elements  and  subgrade  soils." 

Results  from  our  study  also  indicate  that  the  modified  Berggren 
equation  provides  reasonable  estimates  ot  the  trost  depth  when  accurate 
values  of  the  initial  and  boundary  conditions  are  provided.  The  equa- 
tion will  provide  reasonable  estimates  of  frost  penetration  depths  and 
can  be  used  bv  N.IDOT  to  assess  frost  penetration  beneath  various  pavement 
cross  sections  with  differing  thermal  properties. 


Jo 


Table  8.  Initial  and  upper  boundary  conditions 
for  test  s i tes . 


Site 

No 

Air  or  surface 
freezing  index 
(°F-davs) 

n- factor 

Length  of 
season  days 

Ini t l a 1 
t empera  t un 
(*F) 

1 

44b 

0.  33 

44 

15. 

44b 

0.  53 

44 

35. 

3 

44b 

0.55 

44 

35. 

4 

44b 

0.45 

44 

15. 

5 

44b 

0.45 

44 

IS  , 

h 

See  Table  4 

7 

44b 

0.55 

44 

35. 

8 

44b 

0.55 

44 

15. 

g 

See  Table  4 

10 

44b 

0.45 

44 

35. 

11 

See  Table  4 

12 

388 

l .00 

70 

35. 

1 3 

388 

1 . 00 

70 

35. 

14 

388 

1.00 

70 

35. 

IS 

388 

1.00 

70 

35. 

lb 

822 

0.55 

70 

35. 

17 

822 

0.55 

70 

35. 

18 

See  Tab  1 e 4 

19 

822 

0.  55 

70 

35. 

20 

822 

0.45 

70 

35. 

21 

822 

0.45 

70 

35. 

22 

See  Table  4 

23 

822 

0.45 

70 

35. 

24 

822 

0.45 

70 

35. 

25 

822 

0.45 

70 

35. 

2b 

103 

1 . 00 

70 

35  . 

27 

107  3 

0.  55 

70 

35  . 

28 

No  Calculations 

29 

See  Table  4 

30 

107  S 

0.55 

70 

35  . 

NOTF : 


An  n-faotor  of  1.00  Indicates  that  the  surface  troe.'ing  index  was 
used . 
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i i«r  p t ha 
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ma  ik  t mum  ! rout 
wlnt  «*i  . 


^ «i  imtfl  f f*>at  Iri’t  h 

» • h k Difference* 

H. 

) 

M.l 

Id  i>aa 

4.0 

10.0 

3.5 

10.5 

11.5 

10.0 

It 

-5 . 5 

5.0 

4 

.0.0 

1 4 . 

26.0 

10 

u.s 

20.  S 

11.5 

11 

18 . It 

12.0** 

6.0 

l.’ 

18. 0 

12.0 

6.0 

1 1 

18.  S 

12.0 

6 . 5 

14 

16 . 0 

12.0 

4.0 

IS 

ts.o 

11.5 

1.5 

16 

27.  5 

19. 5 

-12.0 

17 

17.0 

17.0 

0.0 

18 

17.0 

44. 0*  * 

-7.0 

19 

18.0 

16.0 

2.0 

20 

40. 5 

12.5 

8.0 

21 

17.5 

13.5 

4.0 

->  ■> 

40.0 

14.0** 

6.0 

21 

40.0 

10.0 

10.0 

24 

11.5 

15.0 

-1.5 

25 

12.0 

35.0 

-3.0 

26 

35.0 

38.5 

-3.5 

27 

.... 

61.0 

- 

28 

*** 

*** 

- 

29 

52.0 

*9.0** 

1.0 

10 

41.5 

55.0 

-11.5 

* 

D 1 f f erence 

= [Measured  - 

Mean  Difference  = 
Standard  Deviation  = 

Calculated  1 . 

1.8 
+ 7.8 

** 

See  T.ihlf 

5 for  calculate 

d frost  depths. 

*** 

Bedrock  .it 

10. 5- in.  depth 

. Frost  penetration  exceeded 

this  depth 

****  Possible  error  In  observed  data. 
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DRAINAGE  LAYERS 

Berg  and  McGaw's  (1978)  results  Indicated  that  the  use  of  stabilized 
or  unstabilized  open  graded  drainage  layers  would  not  significantly 
increase  the  maximum  seasonal  frost  penetration  beneath  pavements  in  New 
Jersey.  To  further  evaluate  effects  of  these  layers  on  maximum  frost 
penetrations,  pavement  profiles  incorporating  drainage  layers  were 
assumed  to  exist  at  four  of  the  test  sites.  Two  of  the  four  sites  used 
Portland  cement  concrete  pavements  and  the  other  two  used  asphaltic 
concrete  pavements.  Both  stabilized  and  unstabilized  drainage  layers 
were  assumed  beneath  the  PCC  pavements  but  only  stabilized  drainage 
layers  were  assumed  beneath  the  asphaltic  concrete  pavements.  Drainage 
layer  properties  were  similar  to  those  used  by  Berg  and  McGaw  (1978). 

All  or  part  of  the  base  course  or  subbase  course  material  at  each  of  the 
sites  was  replaced  by  the  drainage  layer.  The  four  sites  and  the 
properties  of  materials  within  the  hypothetical  pavements  are  shown  in 
Table  10.  Upper  boundary  conditions  and  initial  conditions  are  the  same- 
as  those  shown  in  Tables  5 and  8.  For  Site  No.  22,  Solution  1 conditions 
(Table  5)  were  used  and  for  Site  No.  29,  conditions  of  Solution  2. 

In  Table  11,  calculated  frost  depths  for  the  assumed  pavements 
incorporating  drainage  layers  were  compared  with  computed  frost  depths 
shown  in  Table  9 for  the  pavements  without  drainage  layers.  In  all 
cases,  calculated  frost  depths  for  pavements  with  drainage  layers  were 
equal  to,  or  slightly  less  than,  those  for  similar  pavements  without 
drainage  layers. 


Table  11.  Calculated  frost  depths  for  hypothetical 
pavements  with  drainage  layers. 


Calculated  frost  penetration  (in.) 


Site 

No. 

Without  drainage  layer 

With  drainage  layer 
Stabilized  Unstabilized 

L 

21  .S 

21.5 

20.5 

8 

24.0 

21.5 

- 

22 

34.0 

30.0 

27.5 

29 

49.0 

42.0 

- 

29 


Table  10.  Properties  of  materials  in  hypothetical 
pavements  with  drainage  layers. 


Yd 

K 

s 

C 

L 

Material 

w 

d 

/ Btu  . 

Btu 

, Btu 

property*  (lb/cu  ft) 

(%) 

(ft) 

s~ 

»-h 

rt 

o 

l'u  ft  ®F 

cu  f t 

Site  No . 4 , 

st  ah i 1 l zed 

OGL*  . 

PCC 

140. 

0. 

0.83 

0.  57 

30. 

0. 

OGL*  (stab) 

120. 

4 . 

0.50 

0.51 

25. 

691 . 

G 

125. 

12.8 

0.50 

2.  14 

33. 

2304. 

S 

125. 

14.0 

1.0 

2.05 

35. 

2556. 

Site  No.  4, 

unstab  i 1 Iz 

ed  OGL, 

PCC 

140. 

0. 

0.83 

0.  57 

10. 

0. 

OGI.  (unstab) 

121. 

4.5 

0.  50 

0.40 

26. 

784. 

G 

125. 

12.8 

0.50 

2.14 

33. 

2304. 

S 

125. 

14.0 

1.0 

2.05 

35. 

2556. 

Site  No.  8. 

stab l 1 i zed 

OGL . 

AC 

138. 

0. 

0.79 

0.54 

28. 

0. 

OGL  (stab) 

120. 

4 . 

0.42 

0.51 

25. 

691 . 

5A 

138. 

12.0 

0.42 

3.01 

36. 

2385. 

G-S-M 

130. 

12.5 

1 .0 

2.44 

34. 

2340. 

Site  No.  22, 

stabilized  OGL. 

PCC 

140. 

0. 

0.83 

0.  57 

30. 

0. 

OGL  (stab) 

120. 

4.0 

0.5 

0.51 

25. 

691 . 

IA 

130. 

8.8 

0.5 

1 .98 

31. 

1647. 

S-M 

125. 

10.  5 

1 .0 

1.90 

31 . 

1890. 

S-M 

125. 

9.9 

2.0 

1.88 

29. 

1782. 

Site  No.  22, 

unstablllzed  OGI 

, . 

PCC 

140. 

0. 

0.83 

0.57 

30. 

0. 

OGL  (unstab) 

121. 

4.5 

0.5 

0.40 

26. 

7 

OO 

r 

1A 

130. 

8.8 

0.5 

1.96 

31. 

1647. 

S-M 

125. 

10.5 

1.0 

1.90 

31. 

1890. 

S-M 

125. 

9.9 

2.0 

1.88 

29. 

1782. 

Site  No.  29, 

stabilized  OGL. 

AC 

138. 

0. 

0.25 

0.54 

28. 

0. 

AC 

138. 

0. 

0.33 

0.54 

28. 

0. 

OGL  (stab) 

120 

4.0 

0.50 

0.51 

25. 

691. 

1C 

125. 

8.5 

1.0 

1.67 

29. 

1530. 

G-S-M 

125. 

9.0 

2.0 

1.73 

30. 

1620. 

*Opon  graded 

drainage 

layer . 
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DESIGN  VROST  PENETRATION  DEPTHS 


The  Il>7h-1977  winter  was  ranch  colder  than  normal  in  New  Jersey  as 
well  as  most  of  the  northeastern  part  of  the  P.S.  Maximum  depths  ot 
frost  penetration  were  also  probably  greater  during  the  H7S-1H77  winter 
than  in  more  "normal"  winters.  As  a result  of  this  study  frost  depths 
were  measured  at  several  locations  within  the  state  of  New  Jersey 
during  this  cold  winter;  these  recorded  depths  could  be  used  for  pave- 
ment design  purposes.  Computed  results  indicate  that  frost  depths 
beneath  pavements  containing  open-graded  drainage  layers  will  not  be 
significantly  greater  than  beneath  "conventional"  pavements  in  New 
Jersey. 


Tt  the  designers  choose  to  compute  design  frost  penetration  depths 
for  New  Jersey  pavements,  the  modified  Rerggren  equation  can  ho  used. 
Air  freezing  indexes  from  the  l‘)7b-l<l7  7 winter  could  he  used  with  the 
appropriate  n-factors  reported  herein  (Table  J)  for  the  upper  boundar\ 
condition.  A temperature  of  3S°F  should  be  assumed  for  the  soil  prior 
to  the  onset  of  freezing  temperatures  and  the  length  of  the  freezing 
season  at  a particular  site  should  he  obtained  from  NOAA  (H7P-1977) 
records  nearby.  Thermal  properties  ot  the  soils  and  pavement  materials 
used  herein  can  he  used  in  the  computat  ions . Kersten's  t,  1 qc, c> J values 
for  thermal  conductivity  and  volumetric  latent  heat  of  fusion  are  ade- 
quate for  pavement  design  in  New  Jersey. 


CONCl.l'STONS  AND  RECOMMENDATIONS 

The  modified  Rerggren  equation  can  he  used  to  estimate  frost  pene- 
tration beneath  New  Jersey  pavements.  For  the  most  accurate  estimates, 
values  of  the  initial  and  upper  boundary  conditions  and  the  thermal 
properties  and  moisture  contents  of  the  pavement  materials  must  he 
representative  of  those  at  the  site. 

These  calculations  and  those  hv  Berg  and  McCaw  (10 ’SI  indicate  that 
open  graded  drainage  lovers  will  probably  have  vorv  1 i t t le  effect  on  the 
maximum  seasonal  frost  penetration  beneath  pavements  in  New  Jersey . In 
all  ot  the  cal cu lat ions , however,  it  was  assumed  that  the  pavement 
surface  temperature  over  a drainage  lover  is  the  same  as  that  over  a 
conventional  base  course.  It  was  also  assumed  that  the  moisture  condi- 
tions beneath  a drainage  lover  ore  the  same  as  those  in  a more  conventional 
pavement.  Errors  in  either  or  both  of  these  assumptions  could  cause  a 
somewhat  greater  frost  penetration,  in  the  order  of  10  to  IS*  maximum. 

It  (s  probably  useful  to  consider  the  likelihood  that  these  assume 
tions  are  incorrect . As  discussed  in  the  "Other  Cons iderat ions"  port  ion 
of  the  previous  report  (Rerg  and  McCaw  19781,  the  placement  ot  insulating 
layers  in  pavements  can  cause  icing  condit ions  similar  to  those  oocurt ing 


U 


on  bridge  decks.  A relatively  dry  open-graded  drainage  layer  mav  act 
somewhat  similar  to  an  insulating  layer  and  thereby  cause  a higher  sur- 
face freezing  index.  However,  the  surface  freezing  index  over  a pavement 
incorporating  an  open-graded  drainage  layer  would  have  to  increase  over 
10%  to  cause  a 15%  increase  in  the  computed  frost  depth.  It  is  believed 
that  a change  this  large  is  verv  unlikely. 

Should  the  moisture  conditions  not  be  similar  beneath  "convent ional" 

pavements  and  pavements  containing  an  open-graded  drainage  layer,  both 

the  thermal  conductivity  and  the  volumetric  latent  heat  of  fusion  terms 

in  the  modified  Berggren  equation  would  change.  To  cause  an  increase  of 

15%  in  frost  penetration,  the  ratio  K/L  in  the  equation  would  also  be 

required  to  change  more  than  30%.  Such  a difference  may  occur.  For 

example,  a sandv  soil  having  a density  of  110  lb/cu  ft  and  a moisture 

content  by  dry  weight  of  10%  has  L * 1584  Btu/cu  ft,  equation  (2),  and 

unfrozen  K =1.16  Btu/ft  hr  °F,  equation  (3).  If  the  moisture  content 

for  this  soil  were  15%,  L = 2376  Btu/cu  ft  and  K = 1.29  Btu/ft/hr  °F. 

In  this  case  the  ratio  of  K / L would  varv  bv  about  26%. 

s 

Calculated  frost  depths  beneath  hypothetical  pavements  incorporating 
drainage  layers  are  not  significantly  different  from  those  beneath  more 
conventional  pavements  which  use  dense-graded  base  and/or  subbase  course 
materials.  It  is  recommended,  however,  that  prior  to  widespread  use  of 
open  graded  drainage  layers  in  New  Jersey  highway  pavements  test  sections 
in  at  least  two  different  environmental  conditions  be  evaluated. 

The  structural  adequacy  of  the  open  graded  drainage 
also  be  evaluated. 


lavers  must 
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APPENDIX  A 


Observed  Moisture  Contents  and  Frost  Depths 


Plotted  from  data  obtained  by  the 
New  Jersey  Department  of  Transporta- 
tion during  the  1976-1977  winter. 
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